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Abstract:  This  report  presents  a  research  study  conducted  to  predict  the 
impact  forces  that  occur  when  a  barge  train  impacts  head-on  with  a 
circular  concrete  or  concrete-filled  structure.  These  structures  are  found  at 
the  end  of  lock  approach  walls  as  semicircles  called  bullnoses,  as  circular 
cells,  and  as  mooring  cells.  This  research  was  conducted  using  finite 
element  analysis.  A  detailed  finite  element  model  of  the  bow  of  a  jumbo 
hopper  barge  was  created.  LS-DYNA  was  used  to  perform  impact  simu¬ 
lations  of  the  barge  bow  impacting  differing  diameters  of  cell  structures 
and  at  varying  approach  velocities.  In  each  case  a  “capping  force”  was 
found  caused  by  the  plating  and  internal  structure  acting  as  a  “fuseplug” 
during  the  crushing  of  the  bow. 

Although  research  projects  have  been  performed  focused  on  corner  barge 
impacts  with  lock  walls,  little  research  has  been  conducted  in  the  area  of 
forces  from  head-on  collisions  between  a  barge  and  a  bullnose  or  a  cell.  A 
complementary  R&D  effort  has  been  conducted  by  Dr.  Gary  Consolazio 
and  his  research  associates  at  the  University  of  Florida.  His  research  has 
focused  on  a  barge  impacting  a  bridge  pier.  Comparison  with  their  results 
has  been  made  when  it  was  appropriate. 
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1  Introduction 

1.1  Background:  barge  train  maximum  impact  forces 

Locks  are  a  necessary  structural  feature  found  at  every  dam  within  the  U.S. 
inland  waterways  navigation  system.  This  network  of  rivers  is  an  essential 
component  of  the  nation’s  transportation  infrastructure  system,  a  system 
key  to  national  commerce.  Locks  allow  for  groups  of  barges,  lashed 
together  to  form  barge  trains,  to  negotiate  the  changes  in  river  elevation  at 
the  dams.  One  of  the  loads  applied  to  the  locks  of  the  U.S.  inland  waterway 
system  is  the  head-on  impact  of  a  barge  train  with  the  semicircular  portion 
of  bullnose  structures  or  with  concrete-filled  cell  structures.  (In  this  report 
the  term  bullnose  will  be  used  to  refer  to  either  of  these  two  types  of 
structures  as  the  geometry  of  the  contact  area  is  the  same  for  both.)  The 
photograph  on  the  cover  of  this  report  shows  a  badly  damaged  barge  bow 
after  an  impact  with  the  bullnose  on  the  upstream  lock  approach  wall  at 
Smithland  Lock  and  Dam  on  the  Ohio  River.  These  loads  do  not  happen 
frequently  as  they  are  generally  the  result  of  loss  of  control  of  the  barge 
train.  However,  they  do  occur,  and  the  potential  for  occurrence  always 
exists.  Consequently  this  load  case  is  a  design  consideration  for  bullnose 
and  cell  structures. 

This  research  report  discusses  the  results  of  a  series  of  nonlinear  finite 
element  analyses  computing  the  limiting  impact  force  caused  by  yielding 
and  buckling  of  the  deck  and  skin  plates  and  the  internal  structural  frame 
within  the  center  area  of  the  bow  of  a  jumbo  open-hopper  barge  during  a 
head-on  impact  with  a  bullnose  or  cell  structure.  The  structural  concept 
can  be  explained  as  follows.  Because  of  the  elastoplastic  and  limiting  strain 
material  characteristics  of  steel  combined  with  the  structural  layout  of  the 
deck  and  hull  plate,  the  internal  structural  plates,  and  the  angle  steel  of  the 
internal  trusses,  the  multi-degree-of-freedom  structure  of  the  barge  bow 
provides  a  limiting  force  resistance  when  significant  permanent  structural 
deformation  occurs  during  an  impact  event.  When  placed  in  a  severe- 
impact  environment  in  which  the  bow  begins  to  crush,  the  barge  bow  will 
act  like  a  structural  fuseplug  to  provide  for  a  limiting  impact  force  applied 
to  the  bullnose  or  cell  by  the  barge  train.  Thus,  there  is  an  upper-bound 
force  that  bullnoses  and  cells  will  be  subjected  to  during  a  head-on  impact 
event.  This  report  summarizes  the  research  effort  investigating  the 
magnitude  of  this  limiting  force.  The  barge  used  in  this  study  is  a  rake 
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section,  jumbo  open-hopper  barge  200  ft  long,  35  ft  wide,  and  13  ft  high^ 
at  the  hopper  section  region. 

1.2  Yielding  and  plate  buckling  analyses 

In  the  LS-DYNA  nonlinear  finite  element  analyses  discussed  in  this  report, 
the  impact  areas  at  the  bow  of  the  barge  and  the  approach  bullnose 
(modeled  as  a  nonpenetrating  rigid  structure)  are  brought  into  contact 
with  each  other  at  a  constant  velocity  in  the  numerical  simulation  of  the 
crushing  of  the  impact  area  of  the  barge  bow.  Crushing  of  the  bow  area 
continues  until  36  in.  of  penetration  of  the  approach  wall  into  the  barge 
bow  is  achieved.  Figure  1.1a  shows  an  overhead  view  of  the  barge  bow  at 
initial  contact  with  the  bullnose  in  the  numerical  analysis  using  LS-DYNA; 
Figure  i.ib  shows  an  overhead  view  of  the  barge  bow  after  36  in.  of  pene¬ 
tration  for  one  of  the  numerical  evaluations.  Figure  1.2  shows  the  corre¬ 
sponding  resultant  contact  force  normal  to  the  rigid  bullnose,  as  computed  by 
LS-DYNA  during  the  course  of  this  numerical  analysis.  The  key  observation  is 
that  the  contact  force  between  the  rigid  nonpenetrating  bullnose  and  the  bow 
increases  in  magnitude  with  bullnose  penetration  into  the  barge  up  to  a  peak 
force  value:  the  limiting  contact  force.  Continued  deformation  beyond 
approximately  10  in.  results  (for  Case  1)  in  a  lower  contact  force. 

Figure  1.2  is  an  example  of  typical  resultant  contact  force  normal  to  the 
rigid  (approach)  bullnose  versus  permanent  deformation  computed  at  the 
point  of  first  contact  with  the  barge  in  the  LS-DYNA  nonlinear  numerical 
analyses.  In  this  figure,  the  normal  force  ranges  in  value  from  o  kip  to  a 
maximum  force  of  2,400  kips  at  9.6  in.  of  permanent  normal  deformation 
at  the  impact  point  of  the  bow.  All  LS-DYNA  analyses  were  continued  out 
to  36  in.  of  permanent  displacement  of  the  bullnose  into  the  bow  of  the 
barge. 

1.3  Nonlinear  structural  dynamic  analyses  using  LS-DYNA 

LS-DYNA  (Livermore  Software  Technology  Corporation  2007)  is  a 
general-purpose  transient  finite  element  program  that  is  used  for 
analyzing  complex  structural  dynamics  problems,  such  as  the  crushing  of 
the  impact  area  of  the  bow  of  a  barge  during  a  head-on  impact  with  a 
bullnose  at  the  end  of  a  lock  wall  or  with  a  concrete-filled  cell  structure.  It 
uses  a  central  difference  scheme  to  solve  the  equation  of  motion  in  time. 


1 A  table  of  factors  for  converting  non-SI  units  of  measure  to  SI  units  is  found  on  page  vii. 
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Key  features  used  in  the  analyses  discussed  in  this  report  are  nonlinear 
dynamics,  use  of  a  nonlinear  elastic-plastic  constitutive  model  for  all  A-36 
steel  plates  and  angles  composing  the  barge,  and  use  of  contact  surface 
formulation  for  a  flexible-body-to-rigid-body  contact.  Lagrangian  shell 
elements  are  used  to  model  the  barge  components. 


Figure  1.1.  Overhead  views  of  the  impact  area  of  the  bow. 
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•Case  #1 


Figure  1.2.  Force  versus  displacement  of  the  initial  contact  point  of  the  bow  for  Case  1,  impact 

with  a  20-ft-dia  meter  bull  nose  at  2  ft/sec. 


1.4  Report  contents 

Chapter  2  discusses  the  finite  element  mesh  for  the  bow  of  the  jumbo 
open-hopper  barge  as  well  as  the  assigned  material  properties.  This  mesh 
was  assembled  using  TrueGRID  (Rainsberger  2006)  and  modeled  in  the 
nonlinear  LS-DYNA  analyses.  Material  properties  corresponding  to  A-36 
steel,  obtained  from  tests  conducted  on  plate  steel  from  a  barge,  are 
assigned  in  these  numerical  analyses. 

Chapter  3  reports  on  the  results  of  the  12  LS-DYNA  nonlinear  finite 
element  analyses.  The  diameter  of  the  bullnose  structure,  the  approach 
velocity,  and  the  location  of  first  contact  between  the  barge  and  the 
bullnose  are  varied  among  the  analyses. 


Chapter  4  presents  the  summary,  results,  and  conclusions. 
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2  Finite  Eiement  Mesh  and  Materiai 
Properties 

2.1  Introduction 

This  chapter  discusses  the  development  of  the  finite  element  mesh  and  the 
material  properties  assigned  to  the  structural  steel  composing  the  barge. 
The  object  of  the  nonlinear  finite  element  analysis  using  LS-DYNA  is  to 
compute  the  limiting  impact  force  resulting  from  the  yielding  and  buckling 
of  the  plates  and  internal  structural  framing  of  a  barge  during  a  head-on 
impact  with  a  bullnose  structure;  therefore,  only  the  bow  region  of  the 
jumbo  open-hopper  barge  needs  to  be  modeled  in  the  numerical  analysis 
(Figure  2.1,  raked  section). 


Figure  2.1.  Side  view  of  the  bow  of  a  raked  section,  jumbo  open-hopper  barge. 


The  bow  is  35  ft  wide,  27.7  ft  long,  and  15.5  ft  high.  At  the  intersection  of 
the  inclined  open-hopper  face  with  the  flat  hopper  base,  the  barge  is  13  ft 
high. 
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2.2  Finite  eiement  mesh 

The  mesh  generation  software  TrueGRID  was  used  to  construct  the  finite 
element  mesh  of  the  bow  of  the  jumbo  open-hopper  barge.  Fully 
integrated  LS-DYNA  shell  elements  (*SECTION_SHELL;  ELFORM  i6) 
were  used  in  a  detailed  model  (Figure  2.2)  of  the  jumbo  open-hopper 
barge  bow  structural  steel  deck,  hull,  and  internal  structural  plates,  as  well 
as  the  internal  angle  steel  making  up  the  internal  structural  trusses.  ^  The 
mesh  consisted  of  357,897  nodes  and  353,646  elements.  The  resulting 
mesh  of  the  bow  of  the  barge  is  shown  in  Figures  2. 3-2.5.  Also  identified 
by  color  designation  in  these  figures  are  the  thicknesses  of  the  structural 
deck,  hull,  and  hopper  plates  of  the  bow  of  the  jumbo  open-hopper  barge. 
The  headlog  and  elliptical  impact  corner  sections  of  the  bow  are  shown  in 
yellow  in  Figures  2.3  and  2.5.  The  rounded  hull  plate  connecting  the  front 
and  side  hull  plates  (in  blue)  is  shown  in  brown  in  Figures  2.2  and  2.3. 
Figure  2.4  shows  the  barge  sliced  longitudinally  along  the  center  line  and 
allows  a  view  of  the  internal  truss  system  within  the  bow  of  the  barge  as 
well  as  the  deck  and  hull  plates.  Seven  internal  trusses  (six  “standard” 
longitudinal  and  one  outboard  longitudinal  truss)  are  contained  within 
this  half-bow  structural  section  model.  The  trusses  are  made  out  of  A-36 
steel  angles.  Figures  2.6  and  2.7  provide  an  internal  view  of  the  actual 
structural  trusses  joining  to  the  vertical  C-channel  and  the  pair  of  hori¬ 
zontal  C-channels  that  provide  rigidity  of  the  headlog  at  the  bow  of  the 
barge.  These  figures  also  show  the  corresponding  finite  element  meshes  of 
these  areas.  Figure  2.8  shows  the  internal  structural  members  with  the 
deck  and  hull  plates  removed.  Figure  2.8a  shows  the  members  as  viewed 
from  above  looking  down  at  the  barge.  The  starboard  side  of  the  barge  is 
on  the  left.  Figure  2.8b  shows  the  internal  structural  members  as  viewed 
from  the  port  side  looking  to  starboard. 


1  LS-DYNA  reduced  integrated  shell  elements  {*SECTION_SHELL;  ELFORM  =  2;  NUMBER  OF 
INTEGRATION  POINTS  =  3)  with  hourglass  control  (^HOURGLASS;  Hourglass  Control  Type=4,  Flanagan- 
Belytschko  stiffness  form)  were  tried  to  save  on  execution  time  in  an  earlier  study  (Ebeling  and  Warren 
2008).  Reduced  integration  resulted  in  execution  times  nearly  one-third  that  for  the  fully  integrated 
elements.  However,  the  results  obtained  varied  significantly  if  the  number  of  processors  used  to  run  a 
problem  was  varied.  Therefore,  the  fully  integrated  shell  element  model  is  used.  The  results  obtained 
from  runs  with  varying  numbers  of  processors  are  observed  to  be  more  consistent. 
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Figure  2.2.  Mesh  of  bow-side  view  of  the  huli  and  hopper  plates. 


Figure  2.3.  Front  eliiptical  corner  view  of  the  deck  plates,  headlog,  and 
elliptical  corner  of  the  bow. 
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Figure  2.4.  View  of  the  bow  of  the  barge  at  the  center-line  cut  of  the  bow. 
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Figure  2.5.  Overhead  view  of  the  bow  of  the  barge. 
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a*  Actual  internal  structual  truss  members 


b.  Corresponding  finite  element  model  and  material 
zone  steel  members,  identified  by  steel  thickness 


Figure  2.6.  Internal  views  of  the  structural  trusses  at  the  bow  as  viewed  from  one 
side  of  the  barge  iooking  iateraily  across  the  barge. 


Figure  2.7.  Another  view  of  the  interior  showing  the  truss  joining  the 
C-channel  of  the  bow. 
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a.  Overhead  view  of  internal  members  (elliptical  corner  plates  in  green  on  right  and  left). 


b.  Internal  members  viewed  from  the  port  side  of  the  barge. 


Figure  2.8.  Internal  structural  members  with  the  deck  and  hull  plates  removed. 
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2.3  Material  modeling  and  material  properties 

With  regard  to  material  modeling  used  for  the  structural  steel,  significant 
yielding  was  expected  to  occur  during  crushing  of  the  deck,  hull  plates,  and 
the  internal  structural  members  in  the  heading  area  of  the  bow  along  with 
the  angle  steel  forming  the  truss  members  affected  by  the  bow  impact.  As 
such,  all  shell  elements  in  the  model  of  the  bow  were  specified  with  a 
plastic-multilinear  material  model  matching  the  stress-strain  data  derived 
from  tests  conducted  on  plate  steel  obtained  for  an  actual  barge.  A  large- 
displacement,  large-strain  formulation  was  used  for  the  shell  elements  in 
the  LS-DYNA  impact  analyses.  Therefore,  the  true  (Cauchy)  stress  and 
true  (logarithmic)  strain  data  were  used  (Figure  2.9).  This  material  model 
for  A-36  steel  was  obtained  from  tests  conducted  at  the  University  of 
Florida  on  standard  18-in.  tension  coupon  (Consolazio  et  al.  2002).  The 
steel  model  for  A-36  steel  has  an  initial  yield  of  36  ksi  (2.480  +  08  N/m^), 
a  yield  strain  of  0.0012,  and  a  failure  strain  of  0.2.  In  the  LS-DYNA 
analyses,  the  material  model  was  specified  as  an  effective  stress  versus 
effective  plastic  strain,  LS-DYNA  material  number  24, 
*MAT_PIECEWISE_LINEAR_PLASTICITY.  Other  material  model 
parameters  specified  as  input  are 

•  E  =  30,000  ksi 

•  Poisson’s  ratio  =  0.33 

•  Yield  stress  =  48  ksi 

•  Rupture  stress  =  86.6  ksi  at  25  percent  strain 

•  20  percent  plastic  strain  failure  (i.e.,  effective  plastic  strain  at  failure 
(EPPF)). 

Note  that  the  rupture  strain  of  20  percent  is  specified  in  these  LS-DYNA 
impact  analyses.  The  implications  are  that  when  the  strain  reaches 
20  percent,  the  strain  integration  point  is  removed  from  the  analysis  by 
LS-DYNA. 

The  LS-DYNA  Material  Model  24  allows  for  the  addition  of  strain  rate 
effects  on  the  material  using  the  Cowper-Symonds  model.  However,  these 
simulations  did  not  implement  any  strain  rate  effects  for  the  A-36  steel 
material  because  the  barge  impact  events  are  not  high-speed  impact 
simulations. 
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Figure  2.9.  True  stress  versus  true  strain  for  A-36  structurai  steel  (from  Consolazio 

et  al.  2002). 
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3  Limiting  impact  Force  Computations 

3.1  Introduction 

This  chapter  discusses  the  results  of  a  series  of  12  nonlinear  finite  element 
analyses  using  LS-DYNA  and  the  353,646-element  mesh  of  the  bow  of  a 
jumbo  open-hopper  barge  with  the  assigned  A-36  steel  material  properties 
discussed  in  Chapter  2.  The  object  of  the  nonlinear  finite  element  analysis 
was  to  compute  the  limiting  impact  force  from  yielding  and  buckling  of  the 
plates  and  internal  structural  framing  at  the  bow  of  the  barge  during  a 
head-on  impact  with  a  bullnose  type  structure  such  as  those  found  at  the 
end  of  a  lock  wall  or  with  a  concrete-filled  cell  structure.  The  nonlinear 
deformations  were  concentrated  in  the  front  half  of  the  bow  region. 
Impacts  were  made  in  such  a  manner  that  identical  load  conditions  were 
evaluated  for  first  contact  occurring  at  two  different  bow  locations.  One 
location  was  the  center  of  the  bow  (only  plating  material  backed  by  trans¬ 
verse  C-channel  members  exists  at  the  center)  as  shown  in  Figure  3.1a. 

The  second  location  was  offset  from  the  barge  center  line  and  in  line  with 
the  first  row  of  truss  structures  to  the  port  side  of  center,  13.1875  in.  from 
the  center  line  of  the  barge  as  shown  in  Figure  3.1b. 


Figure  3.1.  Diagrams  showing  the  two  points  of  first  contact. 
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3.2  Nonlinear  impact  analyses 

Lock  approach  walls  are  designed  for  usual,  unusual,  and  extreme  loads. 
Engineer  Technical  Letter  (ETL)  1110-2-563  (Headquarters,  U.S.  Army 
Corps  of  Engineers  (HQUSACE),  2004)  provides  a  summary  of  the  design 
requirements  for  these  three  load  cases.  The  load  cases  are  defined  in 
terms  of  the  annual  probability  of  exceedance.  The  design  information 
contained  in  HQUSACE  (2004)  is  summarized  in  Table  3.1.  This 
information  also  includes  ranges  for  non-site-specific  velocities  (expressed 
in  local  barge  coordinates)  for  the  three  load  cases  to  use  in  preliminary 
analyses  when  site-specific  traffic  study  results  are  not  yet  available. 

The  impacts  that  are  the  focus  of  this  report  are  for  head-on  collisions.  In  a 
head-on  impact  the  approach  angle  is  90  degrees,  and  the  approach  angle 
column  in  Table  3.1  is  of  no  consequence  for  these  analyses. 


Table  3.1.  Three  design  load  condition  categories,  frequency  of  loadings,  and  typical  ranges 
for  non-site-specific  impact  angles  and  approach  velocities  (from  ETL  1110-2-563). 


Load  Condition 

Annual 
Probability  of 
Exceedance 
(Return  Period) 

Performance 

Criteria 

Velocity,  ft/sec 

Approach 

Angle  0,  deg 

Forward  Vx 

Lateral  Vy 

Usual 

>0.1(1-10 

years) 

No  damage 

0.5-2.0 

0.01-0.1 

5-10 

Unusual 

<0.1  but 
>0.00333 
(10-300 
years) 

Reparable 

damage 

3.0-4.0 

0.4-0.5 

10-20 

Extreme 

<0.00333 
(>300  years) 

Noncollapse 

4.0-6.0 

<1.0 

20-35 

The  concrete-filled  cells  and  bullnose  type  structures  chosen  for  these 
analyses  have  diameters  of  20,  35,  and  50  ft.  The  20-ft  diameter  was 
chosen  as  one  of  the  test  parameters  as  this  is  the  approximate  diameter  of 
the  upstream  bullnose  between  the  two  locks  at  Smithland  Locks  and  Dam 
and  is  a  typical  bullnose  size  at  rigid-wall  lock  and  dams.  A  view  of  the 
upstream  area  of  Smithland  Lock  and  Dam  showing  the  upstream  middle 
wall  bullnose  is  shown  in  Figure  3.2a.  Figure  3.2b  shows  a  close-up  view  of 
the  bullnose. 
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b.  Closer  view  of  the  bullnose,  upper  middle  wall. 


Figure  3.2.  The  upstream  middle  wall  bullnoses  at  Smithland  Locks  and  Dam. 
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The  proposed  flexible  approach  walls  at  Lock  and  Dams  22  and  25  consist 
of  precast  concrete  beams  supported  by  concrete-filled  steel  cells  with  a 
diameter  of  35  ft.  The  cells  at  each  end  of  the  guide  wall  have  a  diameter  of 
50  ft.  These  dimensions  are  typical  of  flexible  guard  wall  and  guide  wall 
designs.  The  cells  at  the  ends  of  the  guide  walls  are  the  cells  that  can  be 
impacted  by  head-on  collisions.  Therefore  50-ft-diameter  structures  were 
chosen  as  another  test  parameter.  Figure  3.3  presents  the  proposed  wall 
design  for  Lock  and  Dam  22. 

The  authors  felt  that  analysis  of  a  third  intermediate  structure  size  would 
be  beneflcial.  Analysis  of  more  than  two  diameters  could  reveal  relation¬ 
ships  between  forces  and  diameter  sizes  more  readily  than  an  analysis  of 
only  two  diameters.  A  decision  was  made  to  include  analysis  of  impacts 
with  a  35-ft-diameter  structure  as  well. 

All  impacts  were  made  with  the  center  barge  impacting  the  bullnose  head- 
on  in  the  analyses  summarized  in  this  report.  That  is,  the  direction  of  the 
velocity  vector  was  along  the  longitudinal  line  of  the  barge,  which  made 
flrst  contact  with  the  most  upstream  point  of  the  bullnose  as  shown  in 
Figure  3.4b  or  with  the  cell  as  shown  in  Figure  3.4a.  For  both  of  these 
cases  the  velocity  of  the  barge  train  normal  to  the  structure  was  equal  to  Vx 
(using  the  Table  3.1  notation)  and  had  no  Vy  component. 
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Figure  3.4.  Barge  train  indicating  movement  in  the  x-direction  only  at  impact. 


Using  the  maximum  forward  velocities  from  Table  3.1  for  the  Usual  and 
Extreme  cases,  the  three  diameters  of  the  impacted  structures  and  the  two 
different  impact  locations  (barge  center  line  and  in  line  with  the  truss),  a 
combination  of  12  sets  of  approach  velocities,  structure  diameters,  and 
impact  locations  were  selected  for  analysis.  Table  3.2  lists  the  structure 
diameters,  approach  velocities,  and  impact  locations. 


Table  3.2.  Parameters  that  are  the  basis  of  the  twelve  test  conditions. 


Diameter 

ft 

Approach  Velocity 
ft/sec 

Impact  Location 

20 

2 

Center  of  Bow 

35 

6 

Truss 

50 

3.3  Results 

This  section  summarizes  the  computed  results  for  the  12  nonlinear  finite 
element  analyses.  The  deck,  the  hull  and  internal  plates,  and  the  internal 
structural  members  all  developed  elastic  and  then  plastic  strains  within 
the  steel  members  as  the  rigid  structure  contacted  the  barge  bow  and  then 
displaced  the  contact  area  of  the  bow. 

Figure  1.1  shows  the  overhead  views  of  the  barge  at  initial  contact  with  the 
structure  and  after  36  in.  of  penetration  for  Case  1.  Case  1  is  for  the  usual 
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velocity  of  2  ft/sec  impacting  a  20-ft-diameter  bullnose/cell.  Impact 
location  for  Case  1  was  at  the  bow  center  line  between  trusses.  Figure  1.2 
shows  the  force  normal  to  the  approach  wall  versus  displacement  of  the 
bow’s  initial  point  of  contact  for  this  case.  This  force  represents  both  the 
normal  contact  force  applied  to  the  approach  wall  and  the  force  imparted 
by  the  wall  onto  the  multi-degree-of-freedom  structural  system  of  the 
barge.  The  normal  force  ranges  in  value  from  o  kip  to  a  maximum  force  of 
2400  kips  at  9.6  in.  of  permanent  normal  deformation  at  the  initial  contact 
point.  The  LS-DYNA  analysis  of  Case  1  continued  until  there  was  36  in.  of 
permanent  deformation  of  the  rigid  wall  into  the  impact  area  of  the  bow. 
The  nonlinearity  in  Figure  1.2  results  from  (1)  the  nonlinear  true-stress 
versus  true-strain  relationship  for  A-36  steel  (Figure  2.9)  and  (2)  the 
removal  of  strain  integration  points  and  elements  when  an  effective  plastic 
strain  of  20  percent  is  achieved  within  the  shell  elements  used  to  model 
the  structural  members  and  plates  during  the  course  of  the  crushing 
analysis  of  the  bow. 

Figure  3.5  shows  the  distribution  of  effective  plastic  strain  at  a  permanent 
deformation  of  9.6  in.  at  the  point  of  impact  with  the  bow  for  Case  1.  Note 
the  “crease”  or  standing  wave  of  deformation  in  the  deck  plate  in  the  upper 
right  picture  as  well  as  the  dimpling  and  creasing  of  the  underside  plates 
just  below  the  vertical  bow  plate  at  this  level  of  deformation.  Figure  3.6 
shows  a  close-up  view  of  the  distribution  of  effective  plastic  strain  at  the 
impact  area  (i.e.,  the  bow  center).  Note  the  erosion  (i.e.,  removal)  of 
elements  by  LS-DYNA  at  the  intersection  of  the  vertical  bow  plate  and  the 
bottom  plates;  this  results  from  an  effective  plastic  strain  of  0.2  (i.e., 

20  percent)  being  achieved  by  shell  elements  within  this  zone.  Figure  3.7 
shows  the  distribution  of  von  Mises  stress  at  a  permanent  deformation  of 
9.6  in.  at  the  point  of  impact  for  Case  1. 
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Figure  3.5.  Distribution  of  effective  piastic  strain  (in  decimal  fraction)  at  9.6  in.  per¬ 
manent  deformation  and  force  versus  displacement  plot  of  the  first  contact  point  of 
the  bow  for  Case  1,  usual  loading  condition  at  the  barge  center  line  and  impact  with 
a  20-ft-diameter  bullnose  structure. 


Figure  3.6.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  9.6  in.  of  per¬ 
manent  deformation  and  erosion  of  elements  in  the  impact  area  for  Case  1,  a  usual 
loading  condition  impacting  a  20-ft-diameter  structure. 
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Figure  3.7.  Distribution  of  von  Mises  stress  at  9.6  in.  permanent  deformation  and 
force  versus  dispiacement  plot  of  the  first  contact  point  of  the  bow  for  Case  1,  usual 
loading  condition  at  the  barge  center  line  and  impact  with  a  20-ft-diameter  bullnose 

structure. 

Figure  3.8  shows  the  distribution  of  effective  plastic  strain  at  a  permanent 
deformation  of  36  in.  at  the  point  of  impact  with  the  bow  for  Case  1.  Note 
the  multiple  or  standing  waves  of  deformation  and  the  large  amount  of 
bending  in  the  deck  plate  in  the  upper  right  picture  as  well  as  the  dimpling 
and  creasing  of  the  underside  plates  just  below  the  vertical  bow  plate  at 
this  level  of  deformation.  Figure  3.9  shows  a  close-up  view  of  the  distri¬ 
bution  of  effective  plastic  strain  at  the  impact  area  (i.e.,  the  bow  center). 
Note  the  erosion  (i.e.,  removal)  of  elements  by  LS-DYNA  at  the  inter¬ 
section  of  the  vertical  bow  plate  and  the  bottom  plates  as  a  result  of  an 
effective  plastic  strain  of  0.2  (i.e.,  20  percent)  being  achieved  by  shell 
elements  within  this  zone.  Figure  3.10  shows  the  distribution  of  von  Mises 
stress  at  a  permanent  deformation  of  36  in.  at  the  point  of  impact  for 
Case  1.  Because  of  the  nonlinear  material  response  and  the  erosion  of 
elements,  the  normal  force  at  36  in.  of  permanent  deformation  shown  in 
Figure  1.2  is  approximately  1,400  kips,  a  significant  reduction  from  the 
maximum  force  of  2,400  kips  at  9.6  in.  of  permanent  normal  deformation 
at  the  impact  area  at  the  center  of  the  bow. 
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Figure  3.8.  Distribution  of  effective  piastic  strain  (in  decimai  fraction)  at  36  in.  per¬ 
manent  deformation  and  force  versus  displacement  plot  of  the  first  contact  point  of 
the  bow  for  Case  1,  usual  loading  condition  at  the  barge  center  line  and  impact  with 
a  20-ft-diameter  bullnose  structure. 


Figure  3.9.  Distribution  of  effective  plastic  strain  (in  decimal  fraction)  at  36  in.  of  per¬ 
manent  deformation  and  erosion  of  elements  in  the  impact  area  for  Case  1,  a  usual 
loading  condition  impacting  a  20-ft-diameter  structure. 
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Figure  3.10.  Distribution  of  von  Mises  stress  at  36  in.  permanent  deformation  and  force  versus  dis¬ 
placement  plot  of  the  first  contact  point  of  the  bow  for  Case  1,  usual  loading  condition  at  the  barge 
center  line  and  impact  with  a  20-ft-diameter  bullnose  structure. 


Figures  3.11-3.16  summarize  the  force  normal  to  the  wall  versus 
displacement  for  the  12  LS-DYNA  cases.  The  computations  are  made  for  o- 
36  in.  of  permanent  deformation  normal  to  the  structure  at  the  point  of 
initial  contact.  The  12  cases  in  Table  3.3  are  presented  in  groups  of  two  per 
figure.  Odd-numbered  LS-DYNA  analysis  cases  correspond  to  impacts  at 
2  ft/sec  while  even-numbered  LS-DYNA  analysis  cases  correspond  to 
impacts  at  6  ft/sec.  Each  figure  shares  a  common  structure  diameter  (20, 
35,  or  50  ft)  and  initial  impact  location  on  the  bow  of  the  barge  (center  line 
or  in  line  with  a  truss). 
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•  Case  #1 
■  Case  #2 


Figure  3.11.  Force  versus  displacement  curves  for  Cases  1  and  2. 


The  force  versus  displacement  curves  for  Cases  i  and  2  are  shown  in 
Figure  3.11.  Cases  1  and  2  are  impacts  with  a  20-ft-diameter  bullnose. 

Case  1  impacts  at  2  ft/sec  while  Case  2  impacts  at  6  ft/sec.  For  both  cases 
there  are  two  distinct  peak  values  that  are  notable.  For  Case  1  a  peak  value 
of  1897  kips  occurs  at  1.90  in.  of  displacement.  This  is  assumed  due  to  the 
crushing  of  the  vertical  plating  on  the  front  of  the  headlog.  After  this 
resistance  is  encountered  from  the  structural  truss  members,  the  force 
increases  to  a  maximum  force  of  2414  kips  at  9.65  in.  of  displacement. 
Case  2  displays  an  almost  identical  first  peak  of  1902  kips  at  1.71  in.  of 
displacement  and  a  maximum  value  of  2316  kips  at  7.54  in.  of 
displacement. 
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Figure  3.12.  Force  versus  displacement  curves  for  Cases  3  and  4. 


Cases  3  and  4  are  impacts  with  a  35-ft-diameter  bullnose.  Case  3  impacts 
at  2  ft/sec  while  Case  4  impacts  at  6  ft/sec.  The  force  versus  displacement 
curves  for  these  cases  are  shown  in  Figure  3.12.  Cases  3  and  4  have  very 
similar  initial  peak  values  of  2049  and  2096  kips  at  displacements  of  1.37 
and  1.41  in.,  respectively.  The  initial  peak  values  are  due  to  the  crushing  of 
the  vertical  headlog  plates.  As  the  structural  truss  members  are  crushed, 
the  maximum  force  of  2483  kips  for  Case  3  is  encountered  at  11.62  in.  of 
displacement.  Likewise  Case  4  shows  a  maximum  force  of  2577  kips  at 
11.27  in-  of  displacement  caused  by  the  crushing  of  the  structural  truss 
members. 
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Figure  3.13.  Force  versus  displacement  curves  for  Cases  5  and  6. 


The  force  versus  displacement  curves  for  Cases  5  and  6  are  shown  in 
Figure  3.13.  Cases  5  and  6  are  impacts  with  a  50-ft-diameter  bullnose. 

Case  5  impacts  at  2  ft/sec  while  Case  6  impacts  at  6  ft/sec.  For  Case  5  the 
peak  value  of  2159  kips  occurs  at  1.18  in.  of  displacement.  This  is  assumed 
due  to  the  crushing  of  the  vertical  plating  on  the  front  of  the  headlog.  After 
this  resistance  is  encountered  from  the  structural  truss  members,  the  force 
increases  to  a  maximum  force  of  2682  kips  at  4.80  in.  of  displacement. 
Case  6  displays  an  almost  identical  first  peak  of  2174  kips  at  1.07  in.  of 
displacement  and  a  maximum  value  of  2836  kips  at  9.58  in.  of 
displacement. 
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Figure  3.14.  Force  versus  displacement  curves  for  Cases  7  and  8. 

Cases  7  and  8  are  impacts  with  a  20-ft-diameter  bullnose  structure  with 
the  center  of  the  impact  occurring  in  line  with  a  structural  truss  rather 
than  at  the  center  line  of  the  barge  as  in  Cases  i  through  6.  The  force 
versus  displacement  curves  for  Cases  7  and  8  are  shown  in  Figure  3.14. 
These  curves  do  not  exhibit  a  decrease  in  force  after  the  crushing  of  the 
headlog  plates.  The  force  continues  to  increase  because  of  the  crushing  of 
the  truss  members  in  line  with  the  impact  until  a  peak  force  is  achieved. 
For  Case  7  the  first  peak  force  is  1984  kips  at  a  displacement  of  4.66  in. 
The  maximum  force  due  to  further  structural  crushing  is  2071  kips  at 
14.44  in-  of  displacement.  For  Case  8  the  first  peak  force  is  1836  kips  at 
4.31  in.  of  displacement.  The  maximum  force  from  structural  crushing  is 
2054  kips  at  13.85  in.  of  displacement. 
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•  Case  #9 

•  Case  #10 


Figure  3.15.  Force  versus  displacement  curves  for  Cases  9  and  10. 


Cases  9  and  lo  impact  with  a  35-ft-diameter  bullnose  structure  also  in  line 
with  a  truss  of  the  barge.  The  force  versus  displacement  curves  for  Cases  9 
and  10  are  shown  in  Figure  3.15.  The  crushing  of  the  truss  structure  in  line 
with  the  impact  results  in  a  first  peak  value  for  Case  9  of  2207  kips  at 
3.40  in.  of  displacement.  Consequent  crushing  of  the  structural  members 
yields  a  maximum  force  of  2308  kips  at  9.98  in.  of  displacement  for 
Case  9.  Case  10  results  indicate  a  first  peak  value  of  2194  kips  at  2.94  in.  of 
displacement  and  a  maximum  force  of  2265  kips  at  16.35  in-  of 
displacement. 
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•  Case  #11 

•  Case  #12 


Figure  3.16.  Force  versus  displacement  curves  for  Cases  11  and  12. 


Cases  11  and  12  impact  with  a  50-ft-diameter  bullnose  structure  also  in 
line  with  a  truss  of  the  barge.  The  force  versus  displacement  curves  for 
Cases  11  and  12  are  shown  in  Figure  3.16.  The  crushing  of  the  truss  struc¬ 
ture  in  line  with  the  impact  results  in  a  first  peak  value  for  Case  11  of  2347 
kips  at  2.34  in.  of  displacement.  Not  much  drop-off  in  forces  occurs  after 
this.  The  authors  feel  that  this  is  due  to  the  greater  diameter  bullnose 
structure  engaging  the  additional  structural  truss  members  earlier  than 
the  smaller  diameter  structures.  Consequent  crushing  of  the  structural 
members  yields  a  maximum  force  of  2513  kips  at  7.71  in.  of  displacement 
for  Case  11.  Case  12  results  indicate  a  first  peak  value  of  2399  kips  at  2.20 
in.  of  displacement  and  a  maximum  force  of  2521  kips  at  12.31  in.  of 
displacement. 


The  maximum  values  indicated  are  summarized  and  shown  in  Table  3.3. 
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Table  3.3.  Maximum  values  indicated  during  head-on  barge  impacts. 


Case 

Bull  nose 

Diameter 

Forward  Velocity 
Vx 

Initial  Contact 

Location  on 
Bow 

Maximum 

Force 

kips 

Displacement  at 
Maximum  Force 

in. 

First 

Peak 

Force 

kips 

Displacement  at 
First  Peak  Force 

in. 

in. 

ft 

in. /sec 

ft/sec 

1 

240 

20 

24 

2 

Center 

2414 

9.65 

1897 

1.90 

2 

240 

20 

72 

6 

Center 

2316 

7.54 

1902 

1.71 

3 

420 

35 

24 

2 

Center 

2483 

11.62 

2049 

1.37 

4 

420 

35 

72 

6 

Center 

2577 

11.27 

2096 

1.41 

5 

600 

50 

24 

2 

Center 

2682 

4.80 

2159 

1.18 

6 

600 

50 

72 

6 

Center 

2836 

9.58 

2174 

1.07 

7 

240 

20 

24 

2 

Truss 

2071 

14.44 

1984 

4.66 

8 

240 

20 

72 

6 

Truss 

2054 

13.85 

1836 

4.31 

9 

420 

35 

24 

2 

Truss 

2308 

9.98 

2207 

3.40 

10 

420 

35 

72 

6 

Truss 

2265 

16.35 

2194 

2.94 

11 

600 

50 

24 

2 

Truss 

2513 

7.71 

2347 

2.34 

12 

600 

50 

72 

6 

Truss 

2521 

12.31 

2399 

2.20 

3.4  Comparison  with  an  eariier  finite  eiement  modei 

At  the  time  these  tests  were  performed,  the  authors  decided  to  compare 
the  results  from  one  of  the  test  conditions  from  a  recently  developed 
(2009)  finite  element  model  used  to  predict  forces  encountered  by  a 
jumbo  hopper  barge  impacting  bridge  piers  with  results  that  could  be 
obtained  with  the  barge  bow  model  described  in  this  report.  To  that  end,  a 
test  was  conducted  with  the  following  parameters: 

•  4-ft-diameter  round  impact  structure 

•  0.4  in./sec.  impact  velocity 

•  Impact  at  the  center  of  the  bow. 

In  the  earlier  Consolazio  and  Cowan  (2003)  nonlinear  finite  element 
analyses  of  a  barge  impacting  a  round  (4-ft-diameter)  rigid  object  at 
0.4  in./sec.,  a  peak  force  of  about  715  kips  was  observed  at  12  in.  of 
deformation.  The  LS-DYNA  analysis  provided  a  peak  force  of  1290  kips  at 
5  in.  of  deformation  (Figure  3.17).  The  authors  consider  this  a  significant 
difference.  The  Consolazio  and  Cowan  work  was  performed  using  ADINA 
rather  than  LS-DYNA  for  the  analyses.  There  are  the  obvious  concerns 
regarding  differences  between  the  barge  models  caused  by  the  size  and 
layout  of  some  of  the  structural  angles  and  plates. 
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Figure  3.17.  Force  versus  displacement  curve  for  4-ft-diameter  impactor. 


Also  the  finite  element  model  presented  in  Chapter  2  of  this  report  uses 
continuous  joining  of  the  fully  integrated  shell  elements  where  parts  (i.e., 
angles,  plates,  etc.)  join  together.  Details  on  the  finite  element  model  were 
sparse  in  the  2003  Consolazio  and  Cowan  paper.  It  was  not  known  if  the 
ADINA  model  joined  parts  in  this  manner  or  attempted  to  model  spot 
welds,  which  may  result  in  smaller  impact  forces,  depending  upon  the 
impact  angle  and  the  orientation  of  the  structural  member.  In  reply  to 
inquiries  by  the  authors  (e-mail  message  dated  10  February  2009)  Dr. 
Gary  Consolazio  suspected  that  “the  differences  are  most  likely  due  to  the 
resolution  of  the  2003  finite  element  model  and  to  aspects  of  the  dynamic 
stabilization  procedures  used  in  ADINA.  ...In  our  current  generation  of 
barge  models,  all  members  (angles,  channels,  etc.)  are  either  fully 
discretized  using  shell  elements  or  use  beam  elements  with  integrated 
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cross-sections  and  nonlinear  material  properties.  In  the  ADINA  model 
from  2003,  portions  of  the  finite  element  model  utilized  simpler  resultant 
beam  elements  for  efficiency.  In  addition,  we  now  use  only  fully  integrated 
shell  elements  in  our  models  as  you  do.”  The  most  recent  Consolazio  head- 
on  impact  research  results  for  his  latest  generation  barge  model  are 
summarized  in  Consolazio  et  al.  (2009)  in  which  LS-DYNA  is  used  to  solve 
for  the  nonlinear  force  versus  displacement  of  a  single  jumbo  open  hopper 
barge  impact  with  a  bridge  support.  By  Equation  8  and  by  Figure  8.c  in 
this  Consolazio  et  al.  (2009)  paper,  a  peak  force  equal  to  1125  kips  was 
computed  during  an  impact  with  a  4-ft-diameter  (round)  rigid  structure. 
This  is  much  more  consistent  with  the  LS-DYNA  computed  peak  force 
result  of  1290  kips  in  this  report. 

3.5  Comparison  with  AASHTO  specification  of  impact  force  as  a 
function  of  crush  depth 

In  1991,  the  American  Association  of  State  Highway  and  Transportation 
Officials  (AASHTO)  published  a  report  that  provided  a  guide  specification 
for  waterborne  vessel  collisions  with  highway  bridges  (AASHTO  1991). 
That  report  specified  a  means  to  compute  depth  of  bow  damage  as  a 
function  of  kinetic  energy  of  the  barge  train  (immediately  prior  to  impact) 
and  collision  impact  force  as  a  function  of  the  depth  of  bow  damage.  This 
approach  assumes  the  kinetic  energy  translates  directly  to  permanent  bow 
deformation.  That  specification  was  taken  from  research  conducted  by 
Meir-Dornberg  (1983).  The  AASHTO  equation  to  compute  bow  damage  is 


14 


KE  2 
5672 


10.2 


(3.1) 


where: 

aB  =  barge  bow  damage  depth  in  feet 
KE  =  barge  impact  energy  in  kip-feet 
Rb  =  ratio  of  Bb  /  35 
Bb  =  barge  width  in  feet. 
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The  AASHTO  specification  further  states  the  head-on  collision  force  shall 
be  determined  in  the  following  manner: 


4112. (a,). (i?,) 

if  Qg  <  0.34  ft 

(3.2) 

[1349  +  110. (a,)].(J«.) 

if  Og  >  0.34  ft 

(3.3) 

where  Pb  is  the  static  barge  impact  force  in  feet. 

At  a  meeting  of  the  U.  S.  Army  Corps  of  Engineers  Vessel/Barge  Impact 
Product  Delivery  Team  (PDT)  \  committee  members  requested  that  the 
results  of  this  research  be  compared  with  the  AASHTO  specification.  This 
report  section  is  dedicated  to  this  task.  It  is  important  to  remember  for 
this  discussion  that  the  research  effort  discussed  in  this  report  is  focused 
on  quantifying  the  maximum  force  during  an  impact  caused  by  the 
yielding  of  plates  and  heading  and  the  buckling  of  truss  members  within 
the  bow  region  of  the  lead/impact  barge.  Neither  the  total  mass  of  the 
barge  (and  its  contents)  nor  the  hydrodynamic  effects  of  the  water 
surrounding  the  barge  are  factors  accounted  for  in  this  LS-DYNA 
analytical  model.  The  AASHTO  specification  is  a  means  to  approximate 
the  impact  force  and  depth  of  damage  to  a  barge  based  upon  the  kinetic 
energy  of  the  barge  during  impact.  As  noted  in  Consolazio  and  Cowan 
(2003),  the  AASHTO  specification  computes  the  kinetic  energy  using 
Equation  3.4: 


KE 


29.2 


(3.4) 


where: 


Ch  =  hydrodynamic  mass  coefficient  (factors  the  effect  of  the  water 
surrounding  the  barge) 

W  =  vessel  weight  in  tonnes  (1  tonne  =  2205  lb) 

V  =  impact  velocity  in  ft/ sec. 

Since  the  entire  impact  barge,  its  mass,  and  its  cargo  represented  in  the 
finite  element  model  are  not  available,  the  kinetic  energy  cannot  be 


1  The  meeting  was  conducted  at  the  U.S.  Army  Engineer  Research  and  Development  Center,  Information 
Technology  Laboratory,  in  Vicksburg,  MS,  on  March  11,  2009. 
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computed  in  the  above  manner.  To  rectify  this,  the  potential  energy  for  the 
barge  in  Case  i  of  the  subject  LS-DYNA  analysis  was  determined  by  com¬ 
puting  the  integral  of  the  normal  force  as  a  function  of  the  permanent 
barge  displacement  of  the  bow.  This  yielded  the  potential  energy  possessed 
by  the  crushing  of  the  bow  of  the  barge,  shown  in  Figure  3.18.  The  normal 
force  versus  permanent  deflection  curve  for  Case  1  can  be  seen  in 
Figure  3.11. 
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Figure  3.18.  Barge  bow  potential  energy  for  LS-DYNA  anaiysis  Case  1. 


In  actuality,  during  a  head-on  impact  involving  a  moving  barge  of 
sufficient  velocity  (and  KE)  that  the  plates,  heading,  and  trusses  within  the 
bow  region  of  the  impacting  front  barge  are  crushed,  the  kinetic  energy  of 
the  barge  train  immediately  prior  to  impact  is  assumed  for  this  simplified 
computation  to  be  transformed  into  the  barge  bow’s  potential  energy. 
Using  this  simplifying  assumption,  the  energy  values  shown  in  the  curve  in 
Figure  3.18  from  the  LS-DYNA  Case  1  analysis  were  set  equal  to  the  KE 
term  used  in  Equation  3.1  to  determine  the  equivalent  AASHTO  depth  of 
damage  for  the  impacting  barge  bow,  damage  depth  as,  for  a  head-on 
impact.  Then  the  determined  barge  bow  damage  depth  was  used  to 
compute  the  impact  force  Pb,  according  to  the  AASHTO  specification  using 
Equations  3.2  and  3.3.  The  AASHTO  impact  forces  versus  the  AASHTO 
barge  bow  damage  depth  were  plotted  in  Figure  3.19.  The  LS-DYNA 
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analysis  Case  i  normal  force  versus  displacement  is  also  shown  in 
Figure  3.19  for  comparison. 


—  LS-Dyna  Case  #1  — AASHTO 


Figure  3.19.  Force  versus  displacement  curves  for  LS-DYNA  analysis  Case  1  and  the  corresponding 

AASHTO  specification. 


Four  points  from  the  two  curves  in  Figure  3.19  are  highlighted  for  com¬ 
parison.  The  first  point  selected  is  the  first  peak  force  from  the  LS-DYNA 
Case  1  curve.  The  second  point  is  at  4  in.  displacement  (0.34  ft)  on  the 
AASHTO  curve  (the  point  where  the  slope  of  the  curve  changes).  The  third 
point  corresponds  to  the  point  of  maximum  force  on  the  LS-DYNA  curve. 
And  the  fourth  point  is  the  final  point  on  the  LS-DYNA  Case  1  analysis 
curve,  37.76  in.  The  four  selected  points  and  their  actual  values  are  also 
reported  in  Table  3.4.  The  point  number  column  refers  to  the  particular 
sample  from  the  two  data  sets. 
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Table  3.4.  Comparison  of  LS-DYNA  and  AASHTO  specification  results  for  the  energy 
determined  from  LS-DYNA  Case  1  analyses. 


Point 

LS-DYNA 

Potential  Energy 

AASHTO 

Displacement 

Normal  Force 
kips 

Bb 

Pb 

kips 

in. 

ft 

kip-in. 

kip-ft 

in. 

ft 

51 

1.9018 

0.1585 

1897.4 

2478 

206.5 

2.208 

0.184 

756.8 

85 

3.2003 

0.2667 

1545.4 

4641.6 

386.8 

4.104 

0.342 

1386.6 

254 

9.6501 

0.8042 

2413.8 

16618.8 

1384.9 

14.1 

1.177 

1478.5 

970 

37.7589 

3.1466 

1494.6 

67219.4 

5601.6 

50.2 

4.180 

1808.8 

As  can  be  observed  by  the  study  summarized  using  the  four  points  listed  in 
Table  3.4,  for  an  equivalent  amount  of  impact  energy  the  AASHTO  speci¬ 
fication  always  results  in  a  permanent  deformation  greater  than  that 
obtained  by  the  LS-DYNA  analyses.  Also  the  LS-DYNA  impact  force  for 
this  case  peaks  before  10  in.  of  permanent  deformation  occurs.  Because  of 
the  buckling  and/or  fracturing  of  the  plates  and  internal  structural  mem¬ 
bers,  a  softening  of  the  bow  occurs  and  the  impacting  force  decreases  from 
the  peak  force.  The  LS-DYNA  impact  force  becomes  less  than  the  AASHTO 
specification  at  about  28.5  in.  of  permanent  deformation.  This  is  in  con¬ 
trast  to  the  AASHTO  specification  impact  force  that  is  a  function  of  the 
crush  depth  and  always  increases  as  crush  depth  increases.  Consolazio  et 
al.  (2009)  also  reported  that  peak  impact  forces  occurred  during  their 
LS-DYNA  analysis  of  barge  impacts  with  bridge  support  structures  before 
maximum  crush  depth  was  encountered. 

3.6  Comparison  with  barges  impacted  using  the  14-IVIN  Statnamic 
ioad  device 

ETL 1110-2-563  (HQUSACE  2004)  summarizes  full-scale  crushing 
experiments  conducted  in  New  Orleans,  Louisiana,  at  the  Halter  Gulf 
Repair  facilities  during  21-23  June  2000.  The  experiments  consisted  of 
using  two  jumbo  open-hopper  35-  by  195-ft  barges  that  were  recently 
removed  from  service  on  the  inland  waterways  and  donated  for  the 
experiments.  The  barges  were  impacted  using  the  14-MN  Statnamic  load 
device  owned  by  Applied  Foundation  Testing  of  Green  Cove  Springs, 
Florida.  The  Statnamic  device  is  used  primarily  to  test  the  axial  and  lateral 
capacities  of  piles  and  drilled  shafts.  The  Statnamic  device  used  for  the 
experiments  has  the  capability  to  deliver  up  to  2,400  kips  of  lateral  force  at 
a  time  duration  similar  to  that  of  a  barge  impact.  Figure  3.20  (Figure  F-3 
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in  ETL 1110-2-563)  shows  the  Statnamic  equipment  and  the  experimental 
setup  during  the  crushing  experiments. 


Figure  3.20.  Statnamic  device  and  experimental  setup  during  experiment  (Note:  The  detonation  of  the 
Statnamic  device  had  just  occurred  at  the  time  of  the  photograph)  ((HQUSACE  2004)). 


A  total  of  nine  experiments  were  conducted  on  both  the  barge  corners  and 
headings  (front  face  of  the  barge  above  the  rake)  of  the  two  barges  to 
determine  the  impact  forces  and  deformations  of  the  components.  The 
impact  loads  on  the  barges  ranged  from  400  kips  up  to  1,600  kips  of 
lateral  forces.  Deformations  range  from  no  observable  to  a  foot  of 
displacement.  Figure  3.21  (Figure  F-4  in  FTL 1110-2-563)  shows  the 
crushing  damage  to  the  heading  of  the  barge  under  an  800-kip  force 
applied  between  rake  trusses.  Maximum  deformation  of  the  heading  in 
Figure  3.21  is  approximately  9  in. 
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Figure  3.21.  Crushing  damage  of  barge  headiog  (from  HQUSACE  2004). 


Direct  comparisons  between  the  results  of  the  field  Statnamic  impact 
experiments  and  the  LS-DYNA  analyses  are  not  possible  because 
deflections  at  the  stages  of  loading  by  the  Statnamic  device  were  not  made 
so  as  to  develop  figures  similar  to  3.11  through  3.17.  However,  some 
observations  are  made  by  the  authors  of  this  report: 

•  Table  3.3  peak  force  data  and  the  peak  force  in  Figure  3.17  show  that 
the  peak  force  encountered  decreases  with  the  diameter  or 
equivalently,  the  width  of  the  impactor. 

•  The  Statnamic  impactor  with  a  width  that  is  much  smaller  than  that  of 
a  cellular  structure  and  smaller  than  the  4-ft-diameter  impactor 
discussed  in  Section  3.4  will  yield  smaller  impact  forces. 

•  For  the  Statnamic  impactor,  the  peak  forces  encountered  will  vary 
dependent  upon  whether  the  impact  occurs  centered  on  a  truss  or 
between  trusses.  These  locations  were  not  reported  in  ETL 1110-2-563 
(HQUSACE  2004). 

•  The  lower  force  Statnamic  tests  that  resulted  in  zero  or  nominal 
deformation  likely  did  not  contain  a  sufficient  charge  in  the  Statnamic 
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device  to  develop  an  impact  force  that  can  be  associated  with  peak 
forces  in  the  heading  during  crushing.  This  would  be  analogous  to 
determining  the  force  on  the  “linear”  portion  of  load  versus  deflection 
curve  of  Figures  3.11  through  3.16  for  points  that  are  below  a  point 
identified  as  ‘first  peak  force”  in  Table  3.3. 

•  The  9  in.  of  crushing  damage  reported  in  ETL 1110-2-563  for 
Experiment  9  (Figure  3.21)  may  not  be  associated  with  the  peak 
reported  force  of  800  kips. 

•  The  deflection  (and  corresponding  load)  associated  with  the  first  peak 
force  listed  in  Table  3.3  is  likely  dominated  by  the  crushing  of  the 
heading.  As  deflections  increase,  the  truss  system  is  engaged.  This 
results  in  additional  load-carrying  capacity  and  a  second  peak  force 
value,  identified  as  Maximum  Force  in  Table  3.3.  Without  deflection 
data  (in  the  format  of  a  load  versus  deflection  curve)  for  the  Statnamic 
impact  tests,  it  is  not  possible  to  confidently  determine  with  what  stage 
of  loading  or  level  of  structural  system  damage  the  reported  force  is 
associated. 

•  The  barges  tested  were  at  the  end  of  their  service  life.  The  LS-DYNA 
analyses  were  made  using  geometry  and  material  properties  of  a  new 
jumbo  open-hopper  barge.  Therefore,  higher  impact  forces  are  likely 
for  the  LS-DYNA  analyses. 
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4  Summary,  Results,  and  Conclusions 

4.1  Summary 

This  research  report  discusses  the  results  of  a  series  of  12  nonlinear  finite 
element  analyses  using  LS-DYNA  to  compute  the  limiting  impact  force 
from  the  yielding  and  buckling  of  plates  and  internal  structural  framing  at 
the  impact  point  of  the  bow  of  a  barge  during  its  head-on  impact  with  a 
bullnose  structure  at  the  end  of  a  lock  wall  or  with  a  concrete-filled  cell 
structure.  A  finite  element  mesh  of  the  bow  of  a  jumbo  open-hopper  barge 
consisting  of  357,897  nodes  and  353,646  elements  was  constructed. 
Because  of  the  localization  of  the  nonlinear  response  during  impact  of  the 
bow  of  the  barge  with  the  approach  structure,  only  the  bow  region  of  the 
jumbo  open-hopper  barge  (Figure  2.1)  was  modeled  in  the  numerical 
analyses. 

4.2  Results 

Significant  yielding  occurred  during  crushing  of  the  deck  and  hull  plates 
and  the  internal  structural  members  contained  within  the  impact  zone  of 
the  bow  in  all  12  analyses.  All  12  numerical  analyses  were  carried  out  until 
36  in.  of  permanent  deformation  occurred  at  the  first  contact  point  of  the 
bow  of  the  barge.  As  such,  all  fully  integrated  shell  elements  in  the  model 
of  the  bow  were  specified  with  a  plastic-multilinear  material  model 
matching  the  stress-strain  data  derived  from  tests  conducted  on  A-36  plate 
steel  obtained  for  an  actual  barge.  A  large-displacement,  large-strain 
formulation  was  used  for  the  shell  elements  in  the  LS-DYNA  impact 
analyses  and  the  true  (Cauchy)  stress,  true  (logarithmic)  strain  data 
(Figure  2.9)  specified  for  this  material  model  for  the  A-36  steel.  These  test 
results  were  obtained  from  tests  conducted  at  the  University  of  Florida  on 
a  standard  18-in.  tension  coupon  tested  by  Anderson  and  summarized  in 
Consolazio  et  al.  (2002). 

4.3  Conclusions 

This  technical  report  provides  a  basis  for  establishing  the  “capping  force” 
when  a  barge  train  impacts  a  bullnose  structure  head-on.  The  capping 
force  is  based  on  the  computation  of  the  crushing  force  imparted  during  a 
head-on  impact  of  the  bow  of  a  barge  with  a  bullnose  structure.  When 
placed  in  a  severe  impact  environment  in  which  the  bow  begins  to  crush. 
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the  bow  acts  like  a  structural  fuseplug  and  provides  for  a  limiting  impact 
force  applied  to  the  approach  structure  by  the  barge  train.  The  results  of 
these  12  nonlinear  finite  element  analyses  of  crushing  of  the  barge  bow 
definitively  show  a  limiting  force  that  may  be  transferred  to  a  bullnose  by 
the  barge  train  through  the  front  barge  bow. 

The  magnitude  of  this  limiting  force  ranges  in  value  from  2050  to 
2840  kips  for  the  12  analyses.  The  peak  forces  predicted  were  proportional 
to  the  diameter  of  the  bullnose  that  was  impacted  as  shown  in  Table  4.1.  In 
all  cases  the  peak  force  increased  as  the  diameter  of  the  bullnose 
increased.  This  was  to  be  expected  as  more  resistance  from  the  plating  and 
structural  reinforcement  is  encountered  because  a  larger  cross  section  of 
the  barge  is  engaged  because  of  the  larger  bullnose  diameter. 

Cases  1-6  were  identical  to  Cases  7-12  except  for  the  point  of  first  contact. 
Cases  1-6  made  first  contact  at  the  center  of  the  bow.  Cases  7-12  made  first 
contact  in  line  with  the  truss  just  to  the  port  side  of  the  center  line.  In  each 
instance  the  case  that  impacted  at  the  center  always  resulted  in  a  greater 
peak  force  than  the  corresponding  case  that  made  first  impact  at  truss 
location  (i.e.,  Case  1  yielded  a  larger  peak  force  than  Case  7,  Case  2  was 
greater  than  Case  8,  etc.). 

Therefore,  the  peak  force  predicted  is  dependent  upon  the  size  of  the 
bullnose  impacted  and  the  point  along  the  bow  where  the  impact  occurs. 


Table  4.1  Summary  of  maximum  impact 
forces  during  head-on  impacts. 


Diameter,  ft 

Maximum  Forces  Range,  kips 

Minimum 

Maximum 

20 

2054 

2414 

35 

2265 

2577 

50 

2513 

2836 
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